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ABSTRACT 

The  loss  of  particles  through  a  cusp  of  a  particular 
containment  geometry  utilizing  cusped  magnetic  field  lines 
Is  considered.   A  velocity  space  loss  criterion  analogous 
to  the  loss  cone  In  the  mirror  machine  Is  derived.   The 
effect  of  a  uniform  longitudinal  magnetic  field  perpendicular 
to  the  containing  field  Is  considered  and  a  loss  criterion 
Is  derlvedo   The  effect  of  the  longitudinal  field  Is  to 
decrease  cusp  losses. 
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NYO-9370 

PREFACE 

This  report  presents  some  calculations  done  during 
the  suinmer  of  1956 o   The  results  were  presented  at  the 
Project  Sherwood  meeting  held  in  Berkeley  In  1957 •   A 
short  s\immary  of  the  calculations  appears  in  the  proceed- 
ings of  that  meeting  [l]o 

The  idea  of  using  a  generalized  adlabatlc  invariant 
in  the  calculation  was  suggested  by  Dr.  Harold  Grad  [2,3,k.^, 
I  should  like  to  acknowledge  many  helpful  discussions  with 
Dr,  J.  Berkowitz  and  Dr.  George  Morikawao 
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NYO-9370 

THEORY  OP  CUSPED  GEOMETRIES 
IV,  Particle  Losses  in  Crossed  Fields 

!•   Introduction. 

Consider  the  three-dimensional  cusped  geometry  shown 
in  Figure  1.   In  this  configuration  plasma  is  assumed  to 
be  contained  by  a  magnetic  field  with  field  lines  as  shown 
and  there  is  a  sharp  separation  between  plasma  and  vacuum 
field.   There  exists  an  equilibrium  solution  to  the  hydro- 


FIG.I 


magnetic  free  siirface  problem  and  it  is  well  known  that 
this  geometry  is  hydromagnetically  stable  [3,6,7] • 

In  order  to  consider  particle  losses  we  must  modify 
the  sharp  boundary  between  plasma  and  field  and  assume 
instead  a  bovinding  layer  in  which  plasma  and  field  are  mixed. 
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It  Is  In  this  layer  that  particles  leavlag  the  literlor  are 
turned  around  by  the  field  and  returned  to  the  interior. 
With  this  modification  there  are  then  particle  losses 
through  the  cusps.  The  present  paper  gives  a  report  on  an 
analysis  of  these  cusp  losses. 

We  first  consider  the  case  where  the  interior  region 
is  field  freej  the  particles  move  towards  the  cusps  as 
shown  in  Pigtire  2  and  are  reflected  or  lost  according  to 
some  criterion  to  be  determined.   The  z-axis  is  out  of  the 
plane  of  the 


FIG. 2 


paper.   Let  the  velocity  of  a  particle  be  v  =  (u,v,w)  where 
u,v,w  are  the  x,y,z  components  of  velocity  in  the  coordinate 
system  shown.   We  define 


^   =  tan'  J 


(1) 
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Consider  this  angle  as  the  particle  enters  the  bounding 
layer  at  a  point  where  the  separation  is  d.   We  shall 
determine  a  critical  angle  \|r"(d)  such  that  for  a  particle 
with  \|f(d)  >  V^(d)  it  is  reflected  and  for  \lf(d)  <  ^^{d)    it 
is  lost.  This  velocity  space  loss  region  is  analogous  to 
the  loss 


FIG.  3 


cone  in  the  mirror  machine.  With  such  a  loss  criterion 
we  can  calculate  loss  rates  through  the  cusps. 

We  consider  now  the  case  of  a  B  in  the  interior  region 
only.   Particles  move  in  circxalar  orbits  in  the  x-y  plane, 
i.e,,  they  spiral  about  the  B  lines.   Only  those  particles 
which  are  near  enough  to  the  boundary  to  penetrate  the  layer 
will  move  towards  a  cusp.   There 
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FIG.  4 


is  no  particle  reflection  by  the  converging  field  lines 
until  the  particle  is  close  enough  to  the  cusp  to  penetrate 
both  sides.   It  can  then  be  reflected  or  lost  according  to 
a  criterion  to  be  deterinined. 

The  method  used  in  this  report  to  calculate  losses  is 
to  make  use  of  two  constants  of  the  charged  particle  motion. 
The  first  is  the  speed  which  is  an  absolute  invariant  of 
the  motion;  the  second  is  an  adiabatic  invariant  of  the 
motion.   If  we  consider  the  particle  moticyi  in  the  neighbor- 
hood of  the  cusp  we  see  that  a  particle  penetrating  this 
layer  will  be  turned  around  by  the  magnetic  field  and  re- 
enters the  plasma.   It  then  penetrates  the  layer  on  the 
other  side  and  is  turned  around,  etc.   It  continues  this 
motion  towards  the  cusp  until  it  is  either  reflected  or 
lost.   Now  if  we  consider  one  complete  cycle  of  this  motion 
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and  assume  that  in  the  space  of  one  such  cycle  the  bounding 
surfaces  of  the  plasma  are  parallel,  then  the  motion  is 
periodic  in  velocity  space.   In  configuration  space  it  is 
of  course  periodic  in  one  direction  by  definition  and  is 
iinlformly  displaced  in  the  other  two  directions.   There  is 
a  periodic  motion  in  phase  space  which  can  be  coisidered. 
The  integral  o)  p  .  dq  can  then  be  evaluated  around  this 

closed  curve  in  phase  space,  where  p  and  q  are  the  Hamiltonlan 

-^    -^   e  -^         -» 

variables  of  the  system  i.e.  p  =  mv  +  —  A  ,   where  A  is  the 

c 

magnetic  vector  potentialo   Now  if  we  assume  that  the  separa- 
tion between  the  bounding  surfaces  decreases  gradually  to- 
ward the  cusp,  then  we  assert  that  the  integral  d>p  •  dq  is  an 
invariant  of  the  motion. 

We  shall  use  this  approximate  constant  of  the  motion 
and  the  speed  to  determine  a  reflection  criterion  for  particles 
and  loss  fractions  in  a  manier  analogous  to  the  calculation 
of  the  same  quantities  in  a  magnetic  mirror. 


2o   The  Case  of  no  B   in  the  Plasma  i.e.  B  =  0 

We  consider  the  geometry  as  shown  in  Pigvire  5. 
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FIG. 5 


The   z-axis    is  out  of   the  plane   of  the   paper.      We   ass\iine 
the  bounding   surfaces   to  be   parallel  with   separation  given 
by  d.      In  the   layer  we   assxwie  E  =  0   and  B     =  B  =  constant, 

B     =  B     =  0,      In  the  plasma  E  =  0   and   B  =  0.      The   equations 

y    z 

of  motion  of  a  particle  in  the  layer  are  then 


(2.1) 


mu  =  o 

mv  =  e/c  wB 

mw 


=  -(7c)v  B 


if  we  let 
(2.2) 


mc 
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and  introduce  the  lew  independent  variable 

(2.3)  Q  =  il  t 

then  the  equations  (2„l)  become 

u"  =  o 

(2oI^-)  V'  =  w 

w'  =  -V 

where  the  prime  indicates  d/dO„ 

The  solutions  of  equations  { 2oL|-)  are 


u  =  u 
o 


(2o^)  V  =  w   sin  9  +  V  cos  0 


w  =  =v   sin  9  +  w  cos  9 
o  o 

where  the  initial  conditions  are  u(0)  =u  «  v(o)  =  v  „ 

o  o 

w(o)    =  w   . 

The  particle   trajectories   in   the   layer   are   then  given  by 


X  =    (u  /il  )9  +  X 
o'  o 


(2o6)  y=(v/il)sin9+(w  /-TL   )(l  -    cos   9)    +   y 

"  o  o  o 

z  =   (w  /A  )sin  9   -(v  /rL)(l  =   cos   9)    +   z^ 
o                                o  o 
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where  x(0)  =  x^,  7(0)  =  7^,  z(0)  =  z^o 

In  the  plasma  there  are  no  forces  so  the  particle  moves 
in  a  straight  line  (neglecting  collisions)  with  constant 

velocit7o 

We  now  consider  the  trajector7  shown  in  Figure  5<.   The 
particle  enters  the  la7er  at  the  point  (x^,d/2,z^)  with 
velocit7  components  (u  ,v  ,w  )  at  t  =  9  =  0.   At  the  end 
of  path  1  the  particle  re-enters  the  plasma  at  the  point 
(x,  ,d/2,z-,)  with  velocit7  (i^o>°'^o*%^  where  the  time  of 
path  1  is  given  b7 


^^  =  Q^/r\  =  2/-r\    tan"l    ( -^) 


Hence  we  have 


^1  -  ^o  =  ^o'Tl 


and 


1  -  ^o  =  °  2  ^o/-^ 


At  the  end  of  path  2  the  particles  enters  the  la7er  at 

(Xoo-  4o  Zo)  with  velocit7  (u  o=v^„w^)  and  the  time  of  path 
2*   2    2  -^    o*   o'  o 

2  is 
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^2  "  '^z/^   =  d/v^ 


Then  x^-x,  =  u  T^  and  Zo-z^  =  w  T^,   The  particle  re-enters 
the  plasma  at  the  ead  of  path  3  at  the  point  (x^,-p-,  z^)  with 
velocity  (u  ,v  ,w  ),   The  time  of  path  3  is  given  by 


T_  =  97-n,  =  2/-n.  tan'^(v  /w  ) 
5  5  o   o 


We  have  x-j-Xo  =  u  T,  and  zo-z^  =  2  v  /-^  i- 
5     c.  o   5  5     c.  o 

A  complete  cycle  of  the  motion  is  completed  at  the  end  of 
path  L|.  with  position  (x,  ,  +  ■pjZ,  )  and  velocity  (u  ,v  ,w  ), 
The  time  of  path  1;  is 


T,  =  0,/-n.  =  d/v 


then  X|  -Xo   =  u  T,    and   Zi  -z-v  =  w  T,  .      We    see   that  x,  -x 

u    [2   d/v     +   27i/-^]    since    9.    +   0_    =   2ti,      We    also   have 
o  '    o  '  13 


Zi     -   z     =  w     2d/v 
I4.  o  0^0 


For  this  trajectory  just  considered  we  wish  to  evaluate 


<j)P  • 
We  note  that 


the  integral  ©p  •  dq   around  a  cl'osed  path  in  phase  space. 
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p  =  mv  +  e/c  A 


where  A  is  the  vector  potentialo   We  have  that  B  =  (BpO,0) 
in  the  layer  and  B  =  (0,0,0)  In  the  plasma„   A  is  determined 
by  B  =  curl  A  and  is  satisfied  by  A  =  (0,0,A^)  where  A  (y) 
is  given  by 


B(y  +  d/2)    ;   y  <  -  d/2 


(2.7)        A^  ={ 
z 


0  I  -d/2  <  y  <  d/2 


B(y  -  d/2)    5    y  >  d/2 

We  shall  now  consider  this  motion  in  the  various  phase 
planes*   The  x-motion  is  uniform  with  constant  velocity  u  o 
Consequently,  when  this  constant  motion  is  subtracted  from 
the  x-motlon  the  closed  path  in  phase  space  is  a  point  and 
the  x-component  of  the  iitegral  is  zeroo   The  y-motion  is 
periodic.   The  y-component  of  the  integral  is  then 

<[p  dq   =  imvdy  =  xn/-^  <!   j   v^d9  +  1  v^dQ  +  \    v^d9  +  J  v^d9 


The  z-raotion  is  characterized  by  a  displacement  (zi=z  )  in 
time 


2  d/v  +  2  V-^ 


-  Ik 


so  the  constant  displacement  velocity  Is 

_        w  2d/v 

^2-8)  ^  =  [2d/v  %  27t/^'>-J 

o 

If  this  constant  motion  is  then  subtracted  from  the  z-motion, 
the  closed  path  in  the  phase  plane  is  then  a  circle.   The 
z  component  of  the  integral  is  then  given  by 

<|ip  dq  =  m(})[(w  -  w)+  e/mcA  ](w  -  w)dt 
J  z   Z     J  z 

=  m/-^^!  2j   (w  -  w)^d9  +  I  [w  -  w)4  -^(y  +  d/2)](w-w)de 
(__  o  ® 

+    [(w  -  w  +  -TlCy  -  d/2)](w  -  w)d0( 


Now  let 


u  =  —  [|)p  dq  +  *p  dq  ] 
^   m  J^j   >   J^z  ^z 


Making  use  of  (2o5),  (2.6),  (2.8),  and  Q^,9^,Q^   we  evaluate 
the  above  expression  for  \i.   which  gives 

(2.9)  ti  =  2v  d  +  V~^(vf  +  wf) 

o  o    o 
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This  expression  is  a  generalization  of  the  magnetic 
moinent  for  the  motion  considered,  in  fact  the  second  term 
is  the  magnetic  moment  of  the  particle  in  the  surface 
layer. 

In  order  to  use  the  expressioi  (2»9)  we  consider  d 
to  be  a  gradually  decreasing  function  of  x  with  d(0)  =  0, 
We  assert  that 

(2.10)  ^(d,v,w)  =  constant 

for  a  particle  moving  in  this  region,  where  v  and  w  are 
the  y  and  z  components  of  velocity  of  the  particle  as  it 
enters  the  surface  layer  at  a  point  where  the  separation 
is  d«   We  also  have 


(2.11)  v^  =  u^  +  v^  +  w^  =  constant 


From  (2,10)  and  (2.11)  we  see  that  as  a  particle  moves 
toward  the  cusp,  i.e.  as  d  decreases,  the  x-compoient  of 
velocity,  u,  decreases.   When  the  particle  turns  around  and 
starts  back,  u  =  0,  and  [i   is  given  by 


\i    !„_«  =  2vd  +(x/-n.)v^ 


u=o 
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For  reflection  within  the  device  we  must  have  d  >  0  when 
u  =  Oo   This  condition  is  satisfied  for  particles  with 


(2„i2)  [i  >{%/.^)r 


From  equation  (2.9)  we  have 


|=.6A.((5A)2.^-^-(w/V,2)l/2 

rid 

where  5  =  — rr — ,  the  dimenslonless  variables  expressing  the 
separation.  If  we  let  iJ,"  =(71/-^  HV  ,  the  critical  value  of 
[X  for  reflection,  then  the  above  equation  gives 


(2.13)  (v/V)''"  =  -  5/71  +[(5A)2  +  1  -(w/V)2]l/2 


Now  we  define 


(2. Ill)  ^   =  tan"^  J  =  sln"^  ,  2  ^   2^l/2 

(u  +  V  )  ^ 


\lf  is  the  angle  in  the  x-y  plane  that  the  velocity 
vector  makes  with  the  x-axis.   Then  from  (2.13)  we  have 
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(2,l5)    V"'  =  sin° 


-SA  +  [(5A)^  +  1  -(w/V)^]^^^ 


(1  -(wA)^)^^^ 


Equation  (2.l5)  defines  the  critical  angle  <);"»   At  a  point 
with  separation  d,  particles  entering  the  layer  with  velocity 
(u, v,w)  are  reflected  before  reaching  the  cusp  if  \lf  >  \}/"o 

The  maximum  function  ^"(d)  is  the  one  for  which  w  =  O9 
If  we  consider  only  this  most  conservative  case,  thea  this 
is  the  same  as  considering  the  two  dimensional  problem  in 
the  x-y  planeo   Thus,  the  critical  angle  in  the  two  dimensional 
case  is 


(2.16)    V"  =  sin* 


-5A  +  ((5A)^  +  1)-^^^ 


In  a  report  on  particle  losses  by  Dro  J,  Berkowitz  [5], 
an  expression  for  sin  \lf"  is  derived  by  a  different  methodo 

.;t. 

Both  methods  yield  the  same  expression  for  sin  \|/  „  however,, 

In  PlgTjre  ?_,  sin  \lf"  is  plotted  against  5o   To  compute 
a  loss  fraction,  consider  a  plane  perpendicular  to  the 
X-axis  at  a  point  where  the  separation  is  equal  to  do   If 
we  consider  particles  with  a  fixed  speed,  v,  and  an  isotropic 
distribution  of  velocities  then  the  number  of  particles  passing 
through  this  plane  is  proportional  to  d,  ioeo 

.d/2    Ti/Z 


N--V 


-d/2   \ 


cos\!/d\|;  -~-  V  d 
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Here  we  take  the  angle  \|/  to  be  the  angle  the  velocity 
vector  makes  with  the  x-axis  in  the  x-y  plane  as  the 
particle  passes  through  the  perpendicular  plane.   If  we 
assume  that  particles  passiag  through  this  plane  with 
angle  \|f,  then  enter  the  surface  layer  at  approximately  the 
same  angle  \j;  and  approximately  the  same  separation  d,  then 
the  loss  fraction  of  particles  passing  through  this  plane 
is  given  by 


,d/2    \l/'"'(d) 

r 


)        dy    j       cos\tfd\l» 
^   -d/2       ^o  


^^^^  "  -z/2 — ;i72" 


dy    j       cos\j;d\lf 
-d/2       ^o 


Hence 


(2.17)  L(d)      =    sin  \lf'""(d). 

The  function  sin  \|;"  os  plotted  in  Pig.  7  ^or   the  c^se 
w  =  Oa   For  large  d,  we  have 


(2.18)  L(d)?^1.5  Tr| 


This  gives  an  effective  "hole"  for  losses  of  diameter 
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Instead  of  considering  the  loss  fraction  for  the  case  w  =  0, 
we  can  compute  it  for  the  three  dimensional  case  as  well. 
Let 

u  =  V  cos  \|;  sin  9 
V  =  V  sin  ^  sin  9 
w  =  V  cos  9 


then  equation  (2,15)  becomes 


sin  \|(-"-  = 


sin  9 


'/n  + 


(Vti)^  +  sin 


1/2 


If  we  assume  an  isotropic  distribution  of  velocities  then  the 
loss  fraction  of  particles  passing  through  a  plane  perpendicular 
to  the  X-axis  at  a  point  where  the  separation  is  eaual  to  d  is 
given  by 


,d/2  ,11  p  f\|/-;:-(d,9) 

1  dy  sin     9   d   9      J  cos   ^   d\|; 

L(d)    =    -^/^  '         -^ 


.d/2 


%/2 


]  dy    j       sin   ^0    d9      j  cos   \tf   d\lf 

-d/2  o  o 


Using  the  expression  for  sin  \|f-"-  given  above  we  have 


L(d)  =  2/i:  K   %  +  [(Vti)   +  ij  sln"^ 


(V 


Tl)     +J 


■1/2 


This  expression  for  L(d)  is  slightly  smaller  than  the  result  for 
w  =  0  given  by  equations  (2ol7)  and  (2ol6)<, 
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3,      The  Case  of  a  B   In  the  Plasma. 

We  now  coaslder  the  configuration  of  field  and  plasma 
shown  in  Pig.  6. 


surface  layer 


B. 


T 


plasma 


xB. 


layer 


-•-B. 


■^x 


FIG. 6 


The  component  B  is  perpendicular  to  the  plane  of  the  paper, 
point  out.   The  bounding  surfaces  are  parallel  with  separa- 
tion d.   In  the  layer  we  assume  E  =  0  and  B  =  const., 

B  =  B  =0,   In  the  plasma  E  =  0  and  B  =  const,  =0,  and 
7    z  ^  z 

B  =  B  =0.   The  equations  of  motion  of  a  particle  in  the 

^   y 

plasma  are 


(3.1) 


mu  =  ( ®/c  )  V  B 
z 

mv  =  -(®/c)u  B, 


mw  =  0 
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If  we  let 


eB 


and  Introduce  the  new  independent  variable 

(3.3)        ©  =cJt 

then  the  equations  (3ol)  become 

u'  =  V 

(3.1+) 

v  =  -  u 

w'  =  o 


where  the  prime  denotes  °Tq  » 

The  solutions  of  eauations  (3.ii-)  are 


(3<=5) 


V  =  V  sin  9  +  u  cos  9 
o         o 


V  =  -u  sin  9  +  V  cos 
o         o 


w  =  w 
o 
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where  the  Initial  conditions  are  u(o)  =  u  ,  v(o)  = 


V 


w(o)   =  w   •      The   particle   trajectories   in  the   plasma  are 
given  by 


X  =    (u/tJ)  sin   9  +{v  /uj){l   -    cos    9)    +  x 
o  o  o 

(3.6) 

7  =    (v  A^)  sin  9   -{u  /oJ){l   -    cos   9)    +  y 


z   =    (w  /uJ)    9   +   z^ 
o  o 


where   x(o)   =x    ,    y(o)   =y   ,    z(o)   =   z    a 
o  o'  o 

In  the  plasma  the  particle  has  a  Larmor  radius  given  by 
(3.7)  a  =  iMu^   +  v2)l/2  =  i/^(u^  +  v^)l/2 

The  equations  of  motion  of  a  particle  in  the   surface 
layer  are 


(3.8) 


mu  =  0 


mv  =(^-^c)w  B 


mw  =-(®/c ) V  B 


-  23  - 


If  we  let 


B 

(3.9)  ^  =  B^ 

z 


then  equations  (3o8)  can  be  written 

u'  =  o 
(3.10)         V'  =  bw 
w'  =  -bv 


where  the  prime  denotes  -tqc 

The  solutions  of  equations  (3olO)  are 


u  =  u^ 

^^•■'■■^^         V  =  w  sinb  9  +  V  cosb  9 

o         o 

w  =  -V  sinb  9  +  w  cosb  9 
o  o 

The  particle  trajectories  in  the  layer  are 


X  =  (u  /J)9  +  X 
o'        o 


(3.12)         y  =  (v  /boJ)  sinb  9  +(w  /btj)  (1  -  cosb  9)  +  y 
"  O'  o  o 


z  =  (w  /\>J)    sinb  9  -(v  /bto)(l  -  cosb  9)  +  z 
o  o  o 


Now  consider  a  trajectory  of  the  type  shown  in  Figo  6 
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We  assume  that  we  are  near  enough  to  the  cusp  so  that 


(3.13)         a  >  d  , 

where   a   is  the  Larmor  radius  of  a  particle  when  in  the 

plasma.   Condition  (3.13)  insures  that  the  particle  will 

enter  both  upper  and  lower  surface  layers  for  particles 

with  u  >  0,  i.e.  particles  moving  towards  the  cusp. 

The  particle  enters  the  layer  at  the  point  (x  ,+  ■00  z  ) 

o'      d^      o 

with  velocity  (u  ,v  ,w  ).   At  the  e ^d  of  path  1  the  particle 

000 

enters  the  plasma  at  the  point  (x-,,^,  z^)  with  velocity 

(u  ,  -V  ,w  ),   The  time  of  path  1  is 
0'    o'  o 


Tt  =  9,/c>=  2/hu)ta.n~'^{-v   /w  ), 
X     X  00 


so  we  have 


x,  -  X  =  u  T, 
1    o    o  1 

At  the  end  of  path  2  the  particle  enters  the  lower  layer 
at  the  point  (xp,  -  -^j  Zp)  with  velocity  (UppVp,w  ),  where 


(3.11;) 


Uo  =  u  -  ^  d 

c.  O 


v^  =  -(2(Ju  d  +  v^  -w^d^)-*-/^ 
d.  00 
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The   time   of  path  2   is 


T2   =   62^  ~      /'^si'^    ""'" 


-V  u„   -   u  v„ 
o  2 o   2 

2. ,2 
a    w 


Thus   the   x-displacement   is  given  by 

X2    -   x^   =    -    I/6J   (v^    +    v^) 

At  the  end  of  path  3  the  particle  enters  the  plasma  at  the 
point  (xo,  -  •^j  z^)  with  velocity  {u^,    -   ^2.'   '^o' "      "^^ 
time  of  path  3  is 


Tn  =  ©3/^=  2/bt(J  tan'-'-(-  V2/w^) 


and   the  x-displacement   is  x-,    -  Xp  =   Up   T^.      At   the    end  of 

path  \\  the    cycle    is    completed   at    the   point    (x,  ,    "^  "o*    ziJ 

with  velocity   (u   ,    v    ,    w    ),      The    time   of   path  h,   is 

000 


T^  =   T2,    so 


^Ij.  -  ^3  =   -  ZD  (^o   ■"   ^2^ 
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In  the  trajectory  we  have  considered  the  z  motion  is 
uniform  with  velocity  w  when  the  particle  is  in  the  plasma 
and  differs  from  this  uniform  motion  when  the  particle  is 
in  the  layer.   If  we  make  the  approximation  that  the  z  motion 
is  uniform  with  constant  velocity  w  ,  we  can  subtract  it 
out  completely  when  considering  the  integral 


P  •  dq 


around  a  closed  ciirve  in  phase  space,   Conseauently,  we 
make  this  approximation  and  consider  the  two-dimensional 
problem  in  the  x-y  plane  i.e.  we  set  w  s  o.   With  this 
assumption  we  have 


T  =  T^  =  -K/huJ 


and 


X,  -  X  =  -,  [7:/b{u  +  u„)  -  2(v^  +  v,)] 
h     o    LjJ       O     2        o     d 


The  x-raotion  results  in  a  displacement  x,  -  x^  in  time 
2it/b<j+  29p/^  ,  SO  the  constant  displacement  velocity  is 
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X 


(3.15) 


u  = 


k 


X 


o 


[2n/h     +  29yu;J 


If  this  constant  motion  is  then  subtracted  from  the  x-motion 

the  resulting  path  in  the  phase  plane  will  be  a  closed  curve« 

The  y-motion  is  periodica.   We  now  wish  to  evaluate 

q)p  •  dq  around  the  closed  path  in  phase  space  of  the  trajectory 

consideredo   We  have   p  =  mv  +  e/cA  e   To  determine  A  we  note 

that  B  =  (B  ,  0,  0)  in  the  layer  and  B  =  (0,  0,  B  )  in  the 
X  z 

plasmao   This  magaetic  field  can  be  derived  from  the  vector 
potential 


A 


(A„»  0,  A„) 


where 


(3.16) 


0 


A  ={-B  y 
X   I   z 


0 


y  <  -  d/2 
-d/2  <  y  <  d/2 

J  >   d/2 


and 


(3.17) 


rB^(y  +  d/2)  ; 


^z  = 


0 


V 


B^(y  -  d/2) 


y  <  -  d/2 
d/2  <  y  <  d/2 

y  >  d/2 
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We  then  have,  neglectlag  the  z  motion, 
^  =  1  |)p  .  dq 


©. 


=  ij  ^(u-  u)2  -.  v2]d9  ^^  j  [(u  - 


u]^  +  v^]d© 


+  ^  [  ^/[(u  -  u)  -a;y](u  -  u)  +  v^ 


d9 


9, 

^  j  [(U  -  U)  -Wy](u  -  U)  +  V^l  d9 


Now  making  use  of  (5),  (6),  (7),  (H),  i'^kK    (1^)  ^^^  "the  expressions 

for  9^  ,  9„,  9t,  9,  we  can  evaluate  u,  as  follows: 
1'   2'   3    M-   - 


-,2 


\i  =  VboJ  [vV2   +(Uq    -   u)^]  +  (Uq   -    d/2aJ-u)[(x2- x^)-u  Q^/uj] 
+  ^[  l/SaV©^   +   l/2(v^   -    u^)    sinepCos9p   +  u  v    sin^9p] 


o  o  2  2  o    o 


-^2 


%/huj[v'^/Z  +(U2    -   u)'^]    +(u2   +  d/2(^'   u)[(x^  -   x^)-u  Q^/uj] 


2  2- 


+  l/uj[l/2a  oi'^Q     +   l/2(v2   -   u^)    sin9,  cos9,    +  U2V2sln  9:  ] 

=  n/hiv^/uJ  +  u^d]+l/j  u^v^   -(u^   -6(;d)[v^  +^d(2u^   -^d)]^^^J 

?  1  /?  1 

^r-v    (u„    -wd)    +u^[v^   -■(0H(2u„    -o^n^/ 
2 ,     •       -1         o      2 
a  tOiin  


o      o 


2  ,2 
a  o^ 
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The  above  function  jj,  (ds,u,v)  Is  the  reauired  adiabatic 
Invariant  for  this  type  of  particle  motion.   Again  we 
assume  that  d  is  a  slowly  decreasing  function  of  x  with 
d(0)  =  0,   We  assert  that 

\x    (d,u,v)  =  coasts 


2  2    2    2 
We  have  that  a  a;  =  u  +  v  =  const  for  particle  motion 

in  the  plasma,  and  since  we  are  considering  the  case 

2    2     2 
w  =  0  we  can  assume  V  =  u  +  v  =  const  for  the  eat ire 

2  2    2 
motion.   We  can  then  set  a.  u)     =  V  in  the  expression  for 

[i,   where  V  is  the  speed  of  the  particle. 

When  a  particle  which  is  moving  toward  the  cusp 

turns  around  and  starts  back  u  =  0  and  \i.   is  given  by 

IT]     =  Vbu;[v2]  +  d(v2.wV)^/2  ^  :^  ^„^-l  ujd 
u=o 


For  the  above  expression  to  be  meaningful  we  need  V>dw, 
but  this  is  equivalent  to  (3.13)  which  is  assumed.   Prom 
the  above  we  see  that  the  critical  value  of  [i.   for  re- 
flection within  the  device  is 
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IJ,"  =  Ti/b«i,i  V' 


i.Co  particles  for  which.  \i->[i"   will  be  reflected.   Dividing 
the  equation  for  \i   by  \i"   we  have 


(1)^   I  ^  t<;d 
W^    V  V 


•ji 


bluv   /U    d^    ffV^Z  d/2u    di^f 

r  I/2I 

v,u  (Jd^     .    Ur/VN2   ^d/2u   t<;dv-,  ^ 


Let /-  =  TT—  =  —  and  note  that  sin  ijr  =  -^  and  cos  \|;  =  tt,  then  if 
we  let  \x   =  [i"    the  above  equation  gives  the  equation  which 
defines  the  critical  angle  ^'\ 


I  =  sin  \|;"  +  <(cos\|;"  +  b/71:  sin\|;"  cos\)f' 


(3ol8) 


- '■bA(cos\l/'''"  -  /  )  [sin^\l/^  +/(2  cos\|/^  -/)] 


1/2 


+  ^  sin"-'-  [-sin\l;''''(cos\|f''''-  I)    +   cos\l/'''"[  sin^\|f'""  + /( 2  cos\|;'"'-  L)]-^^^ 


Equation  (3.18)  defines  the  critical  angle  \|/"(d).   At  a  point 
on  the  surface  layer  with  separation  d,  particles  entering 
the  layer  at  an  angle  \|;  are  reflected  before  reaching  the 
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cusp  If  \jf  >  \l;"(d).   Equation  (3»l8)  can  be  solved  numerically 
for  sin  \|/"  as  a  function  of  X  f or  various  values  of  b.   In 
order  to  compare  these  curves  with  each  other  and  with  the 
case  of  no  B   component  of  the  field,  it  is  useful  to  think 
in  terms  of  a  fixed  containiig  field  B   and  plot  the  curves 
against  the  variable  5  =  — rj—  «.  We  note  that  -'^  =  huJ  so 


(3.19)  5=^ 


In  figure  7»  sin  \|f  ,  for  all  the  cases,  is  plotted  against 

5,   In  this  plot  we  cai  see  the  real  effect  of  the  component 

B  ,   In  all  cases,  near  the  cusp  (in  the  neighborhood  of 
z 

d  =  0)  the  curves  approach  the  case  of  no  B  o   Farther  away 

from  the  cusp  the  critical  angle  \|/"  becomes  less,  particularly 

for  the  strong  B   (b  =  1),   In  the  case  of  1^«>  B   field 

z  z 

(b  =  100)  we  see  that  for  sufficient  distance  from  the  cusp 
there  is  even  ai  improvement. 

In  order  to  calculate  the  fraction  lost  of  particles 
passing  through  a  plane  perpendicular  to  the  x-axis,  for  an 
arbitrary  poiat  on  the  x-axis  (i.e.  an  arbitrary  separation 
d),  we  must  also  consider  motion  in  the  region  where  the 
particle  enters  only  one  side  of  the  bounding  layer  or 
neither  side. 
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ro  OJ         — 
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A  particle  which  does  not  enter  either  of  the  boxinding 
layers  makes  a  circular  orbit  in  the  x-y  plane;  consequently, 
it  does  not  move  towards  a  cusp  and  cannot  be  lost. 

Particles  that  enter  only  one  side  of  the  layer  move 
toward  or  away  from  the  cusp  as  shown  in  Pig\ire  I4.,   The 
maximum  distance  that  the  particle  moves  away  from  the  layer 
boundary  on  each  cycle  remains  approximately  the  same  xintil 
the  particle  hits  both  sides.   If  we  assume  the  boundaries 
are  converging  straight  lines  the  motion  is  exactly  periodic. 
Particles  of  one  sign  (charge)  can  only  move  toward  the  cusp 
along  one  side.   For  example,  in  Figure  i|,  is  the  direction 
of  B  is  out  of  the  paper,  positively  charged  psirticles  move 
towards  the  cusp  at  the  right  along  the  lower  side  and  away 
from  the  cusp  along  the  upper  side.   Consequently,  to  find 
the  number  of  particles  lost  in  the  cusp  we  need  only  to 
consider  particles  moving  along  one  side. 

Consider  a  particle  making  this  periodic  motion  towards 
the  cusp  along  the  lower  boundary.  Let  a  particle  enter  the 
plasma  from  the  layer  at  an  angle  \)r  , 


plasma 
field" 


FIG.8 
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The  distance  of  the  particle  from  the  bounding  line  is 
denoted  by  h,  and  h  is  given  by 


(3.20) 


h  =   a[cos\lr    (cos   ©   -   1)    +    sin\lf      sin  9] 


=    a[cos(\|f      -   9)    -   cos\|/    ] 


where  9  =  Cu't  is  measiored  from  when  the  particle  enters  the 
plasma.   The  maximum  distance  of  the  particle  is  H,  given 

by 


(3.21)        H  =  a[l  -  cos^  ] 

o 


A  particle  moving  in  this  manner  will  start  hitting  both 

sides  of  the  bounding  layer  when  the  separation,  d,  decreases 

such  that  d  =  H,   Prom  \j/"(d),  determined  by  equation  (3.l8), 

we  can  evaluate  \|/"(H),  and  consequently,  ^V'  (\|/  )o   ^"(t  )  is 

independent  of  a,  since  we  know  \|f"(  — )  and  thus  \|/"(  — )«   In 

a  a 

Figure  9  we  have  sin  ^"  vs ,  sin  \|;  for  three  values  of  b. 
Prom  the  relation  ^"{"i/  )  we  can  find  a  value  \Jf"  such  that 
for 


\l/  <   \1;         then      \lr   <   \|/ 


and  \1/   >   il/  "      then     \1;   >   \|/ 

Tq  Yg  Yq  Y 
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The  value  sin  \1/"  is  shown  in  Figure  9  as  the  intersection 

between  sin  \\f"    (sin  \l;  )  and  the  line  sin  \|;"  =  sin  \lr  o 

\jf"  is  the  critical  angle  for  particles  movi  ig  toward 

the  cusp  along  one  side  of  the  boundary,  ioe„  for  a 

particle  enteriig  or  leaving  the  layer  at  an  acute  angle  \|/, 

it  will  be  reflected  if  \1/  >  V"  or  lost  if  \}/  <  \|;"o   It  is 

c  c 

evident  that  \|r''"  is  independent  of  d,  since  the  particle 
motion  is  periodic  along  one  sideo 

We  wish  to  calculate  the  fraction  of  the  particles 
passing  through  a  plane  perpendicular  to  the  X-axis  that 
are  lost  through  the  cuspe   We  consider  particles  of  a 
fixed  speed  V,  aid  thus  a  fixed  Larmor  radius,  ao   The 
perpendicular  plane  is  located  as  a  point  where  the  separa- 
tion is  d.   We  assiome  that  we  are  sufficiently  far  from  the 
cusp  so  that  d>2  a,  i.e.  the  particles  can  only  hit  one  side 
of  the  boundary  layer.   Then  the  particles  that  are  lost  are 
only  those  moving  towards  the  cusp  along  one  side„ 

The  iTiraber  of  particles  passing  through  a  plane,  assuming 
a  fixed  speed  and  an  isotropic  distribution  of  velocities 
is  proportional  to 

d/2   %/2 
V^    dy  j  cos\j/d\(f  =  V^  d 

L/2     Jo 
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In  this  integral,  'i/    is  the  angle  that  a  particle  makes  with 
the  X-axis  at  a  point  in  the  plane  with  ordinate  y.   To  find 
the  number  of  particles  lost  we  need  to  calculate  the  integral 


V   I  I  cos  \1/  d\|;  dy 


D 


where  the  domain  of  integration  D  is  priven  by  our  established 
loss  criteriono   We  recall  that  in  the  loss  criterion  the 
angle  \jf  considered  is  the  angle  that  the  particle  path  makes 
with  the  layer  as  it  enters  or  leaves  the  layer,  while  in 
the  above  integral  ijr  is  the  angle  the  particle  path  makes 
with  the  X-axis  passing  through  a  plane  at  ordinate  yo 
Consequently,  we  must  relate  these  angles  to  apply  the 
criterion.   Let  \|/  be  the  angle  at  the  surface  and  let 
(Figure  8)  \lf  be  the  angle  at  a  point  when  the  distance  to 
the  lower  side  is  h.   We  have  from  (3.20) 


h  =   a[cos(\lf      -    9)    -    cos\|/    ] 
o  o 


and  from  \|/  =  tan   —  and 
^         u 


u  =  V  sin  0  +  u  cos  9 
o         o 
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a 


o 
to 
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the  relation 


cos  t  =  sin  \1;   sin  0  +  cos  \lf   cos  9 


=  cos  (\1;^  -  9)  J 


consequently,  we  have 


(3.22) 


h  =   a[  cos   \l/   -    cos   \lf    ] 

o 


Furthermore  If  \1/  =   ^" ,    the  critical  angle,  then 


(3.25) 


h"   =    a[cos   \|f    -    cos    \lf"] 


defines  the  domain  of  integration,  D,   The  domains  are  shown 
In  Figure  10,  for  the  three  values  of  bo   We  then  have 


11 


cos    \|f   d\lr   dy 


t. 


dy 


COS   \|/   dxV  =        h    (\|/)    cos   \|/   d\|r 
-o  o 


A 


=    a 


[cos      \|;   -    cos    \}/'"    COS   \l/ld\|/   =   a/2[\!/"-    sin   ^"    cos   \|/"] 


The   fraction  lost    is   then  given  by 


(3.21;)  L(d)    =   1/2   a/d[^'"'  -    sin   t'"'   cos   \|/'"'] 
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V  _  bV 
then 


To  compare  this  with  the  case  of  no  B  we  note  that  ^   -  ^  -   "77  « 


(3.25)    L(d)  =  V2(^VJ^  [^'^"(b)  -  sin  ^'^{b)    cos  f;'(b)] 


If  we  take  the  values  of  \|;   given  in  Figure  9  we  have 

c 

approximately 


L(d)  =  .18  ^  J  for  b  =  1 

L(d)  =  .52  ^  ^  for  b  =  10 

L(d)  =  '90  ^  I  for  b  =  100, 


which  gives  an  effective  hole   for  particle  losses  of  diameter 


d^  =  .18  ^         ,      b  =  1 


di-  =  .52  —  ,      b  =  10 

^     n. 


d^  =  .90  ^  ,      b  =  100 


This  is  to  be  compared  with 

for  the  case  of  no  B   field  (b  =  oo).   We  see  that  within  the 

z 

framework  of  this  idealized  model  the  addition  of  a  crossed 
field  produces  a  marked  improvement. 

-  39  - 


References 

1.  J.  Killeen,  "Particle  Losses  in  Crossed  Fields",  TID  7536, 
p.  200,  1957  also  UCRL  1|839» 

2.  H.  Grad,  "General  Adlabatic  Theory",  TID  7^20,  p,  1^.02, 
1956,  Gatlinburg  Sherwood  Conference . 

3.  J,  Berkowltz,  K.  0„  Priedrlchs,  H,  Goertzel,  H.  Grad, 
J.  Killeen,  and  E.  Rubin,  "Cusped  Geometries",  Proco  of 
2nd  U.  N,  International  Conference  on  the  Peaceful  Uses 
of  Atomic  Energy,  Vol.  31,  p.  171,  1958. 

Li.,   H.  Grad,  "Theory  of  Cusped  Geometries,  I  General  Survey", 

NYO  7969,  1957o 
5«   J«  Berkowltz,  "Theory  of  Cusped  Geometries,  II  Particle 

losses",  NYO  2536,  1959. 

6.  J.  Berkowltz  et  al,  "Theory  of  Cusped  Geometries, 
III  Computation  of  Cusped  Surfaces",  to  appear, 

7.  H.  Grad,  "Motion  and  Stability  of  Free  Surfaces",  Wash  289, 
p.  Il5,  1955,  Livermore  Sherwood  Conference, 


-1^0  - 


NEW  YORK  UNiytRSITY 

INSTITUTE  OF  MATHEMATICAL  SCIENCES 

LIBRARY 

p5  y/«verly  Place,  New  loi*  3,  N.  Y. 


FEB  23 

1961 

Date  Due 

FtB  1  6  63 

1 

m 

PRINTED 

IN  U.  S.  A. 

,   NYU 

NYO-9370 


c.l 


Killeen 


Theory  of  cusped    ffeoemtries 
IV:    Particle   losses... 


NYU 
NYO-9370 


Killeen 


The o r y  of   cusped   geometrle s. 


TViPflrhinlft    It 


DATE    DUE 


:j^ii:^^^^^E 


BORROWERS    NAME 


N.  Y.  U.  Institute  of 
Mathematical  Sciences 

25  Waverly  Place 
New  York  3,  N.  Y. 


m 

iWili 


